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>• ' Abstract 

gy 

A Coulomb sum rule is derived for the response of nuclei to (e, e') scattering with 

large three- momentum transfers. Unlike the nonrelativistic formulation, the relativistic 

O ! 

Coulomb sum is restricted to spacelike four-momenta for the most direct connection 

with experiments; an immediate consequence is that excitations involving antinucleons, 
e.g., NN pair production, are approximately eliminated from the sum rule. Relativis- 
tic recoil and Fermi motion of target nucleons are correctly incorporated. The sum 
rule decomposes into one- and two-body parts, with correlation information in the sec- 
ond. The one-body part requires information on the nucleon momentum distribution 
function, which is incorporated by a moment expansion method. The sum rule given 
through the second moment (RCSR-II) is tested in 
to be sufficiently accurate for applications to data. 
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through the second moment (RCSR-II) is tested in the Fermi gas model, and is shown 



1 Introduction 

Electrons have proven to be useful probes for the study of nuclear structure, primarily 
because the eN interaction is reasonably well-known and because the electromagnetic in- 
teraction is weak compared to the strong forces which dominate nuclear structure. The 
longitudinal contribution to the nuclear response is of particular interest because of its rela- 
tionship to two-body correlation functions through what is commonly known as the Coulomb 
sum rule (CSR).|IJ With the assumption that the nucleus may be treated as a nonrelativis- 
tic system of interacting nucleons, the extraction of nucleon-nucleon correlations is formally 
direct: this is commonly referred to as the nonrelativistic Coulomb sum rule (NRCSR). In 
practice, however, one finds that to obtain useful information one must extend the exper- 
iments and analysis to energies and momentum transfers sufficiently large that the usual 
nonrelativistic assumptions fail. This will be particularly the case for experiments at the 
Continuous Electron Beam Accelerator Facility (CEBAF) with -Ebeam — 4 GeV. Therefore, it 
has become desirable to have a fully-relativistic formalism both for analyzing experimental 
data and for testing theoretical models of nuclear structure. In this paper, we present an 
extension of the NRCSR which successfully handles certain problems which arise in electron 
scattering from nuclei at large three-momentum transfers. 

Before discussing further the goals of this paper, it will be useful to review the form of 
the NRCSR, as well as its derivation. We begin with the first-order Born (or one-photon 
exchange) approximation for the scattering of ultrarelativistic electrons (|k| » m e ) from 
nuclear targets, in which one can factor the (e, e') differential cross-section into leptonic and 
nuclear parts. After performing a standard separation of this nuclear response function into 
longitudinal and transverse (virtual) photon contributions, the differential cross-section in 
the laboratory frame can be written 
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describes the elastic scattering of ultrarelativistic electrons from a fixed point target with 
charge e and no spin. Here q^= (u>, q) is the four-momentum transfer to the target, 6 is the 
electron scattering angle, E is the electron beam energy and Q 2 = — q 2 = q 2 — uo 2 . The first 
term in the square brackets in ( |1.1| ) is the longitudinal contribution, and is usually expressed 



in terms of the Coulomb response function: 

Wb(w,q) = £ i W(<l)\i}\ 2 S(u -E, + Ek), (1.3) 

/ 

where |i) and |/) denote initial^ and final nuclear states, respectively, and p(q) is the Fourier 
transform of the nuclear charge density operator. 

In general, the operator p(q) should include contributions from both the nucleons and 
(virtual) charged mesons in the target, but it has been conventional in nonrelativistic ap- 
proximations to ignore the latter. Then p(q) is given by the spatial Fourier transform of the 
local nucleon charge density operator: 

p(q) = G E , p (Q 2 )Jd 3 x e^ x ft{x)QJ,(x), (1.4) 

where ^(x) is the local (nonrelativistic) nucleon field-operator, Q = (l-\-Tz)/2 is the charge 
operator in isospin space which projects out protons, and GE, P (q 2 ) is the proton charge form 
factor. (In a purely nonrelativistic treatment of the target, Ge, p (Q 2 ) would be replaced by 
the Fourier transform of the proton charge density _F p (q 2 ).) Charge effects of neutrons are 
small and are usually neglected. 

To proceed to a sum rule, we first define the nonrelativistic Coulomb sum function for 
inelastic scattering at fixed three-momentum transfer q: 
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d 3 xd 3 x> e^ x - x ') E(^ t (x / )Q t ^(x / )l/)(/l^ t (x)Q^(x)|z), (1.5) 
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where the integration is over all energies above the elastic peak at uj = uo e i. Expression (|1.5| ) 
then becomes a sum rule by subtracting the elastic term f—i explicitly and using the closure 
relation 

El/)</l = 1 > (!-6) 

/ 



on the resulting sum over nuclear states | /) . Equation (|1.5|) is then in the form of a Fourier 



transform of the expectation value of four nucleon field operators. To separate out the two- 



1 For notational simplicity, we assume a nondegenerate ground state \i). The results can easily be gener- 
alized to unpolarized targets with J^O. 



body correlation part, we put the operators in normal order (creation operators to the left) 
by anticommutation, using 

Mx),^ (x ' )} = (5(x _ x '). (1.7) 

The sum rule then becomes 

SWq) = Jd 3 x(i\^(x)Q?P(x.)\t) 

+ fd 3 xd 3 x' e iq - (x - x ' ) |(i|^ t (x)Q^ t (x')V'(x / )W(x)|i) 

-(i\$\x)Q$(x)\i)(i\$\^)Q$(x%)\. (1.8) 



The first term is simply the total nuclear charge Z, and the second term can be written in 
terms of a correlation function: 

CW(q) = J d 3 xd 3 x' e^-^^^xW^x'^WQij^li) 

1 lA {i\i,\x)QJ(x)\i}(i\j\x'W(x%}), (1.9) 



Z 
and the elastic form factor of the nuclear target: 



F el ( q ) = G E , p (Q 2 el ) J d 3 x e^ x (i|^(x)Q^(x)|i), (1.10) 

where Ql t = c[ 2 — uj 2 t . Combining (|1.8| )- (|1.10| ), the NRCSR takes the compact form 



S NR (q) = Z + CW(q) - \ |*. (1-11) 

The function Cjv.fi(q) is simply the Fourier transform of the spatial two-proton correlation 
function. It has been defined so that Cnr{<1) = for a system of Z uncorrelated protons. 
Since F e i(q) is measured by elastic (e, e') scattering, the NRCSR ( |1.11| ) provides a direct 
measure of proton-proton correlations. Since both F e i(q) and Cnr{q) are expected to vanish 
in the limit of large |q|, the approach of SVfl(q) ~* Z with increasing |q| is a test of the 
adequacy of the theoretical assumptions. 

Experimental studies of (e, e') carried out in the momentum transfer region up to |q| = 
550 MeV/c over the last decade have left serious questions about the status of the NRCSR. 



For light targets: A = 2 — 4, the measured Coulomb sum Snr{<1) — > Z for the largest mo- 
mentum values reached. || For heavier targets, however, Snr < Z by 20-30% compared to 
theoretical predictions based on independent particle models, [3-8] with larger suppressions 
occurring for larger nuclei. This is despite the fact that Pauli correlations are expected to 
vanish for |q| >2p^~550 MeV/c. We shall not discuss all the possible effects which could be 
involved here, and which have been reviewed elsewhere. |9|, [10] However, one problem which 
is always present in the application of the NRCSR to data is whether the sum defined in 
( |1 . 5|) is in fact saturated by the experimental data: not only is there an absolute relativistic 
limit uj< |q| due to the spacelike nature of virtual photons in electron scattering, but other 
experimental difficulties prevent one from reaching even that limit in practice. We expect 
that extending experiments to higher energies and momenta, such as will be done at CEBAF, 
may improve this situation. 

The main problem dealt with in this paper is the extension of the Coulomb sum rule 
to include the relativity of nucleons in inelastic scattering. This becomes unavoidable for 
momentum transfers |q| ~M, the nucleon mass, even if the nuclear target is essentially non- 
relativistic. This follows because even if the nucleons in a typical nucleus move at essentially 
nonrelativistic velocities, they become highly relativistic after the absorption of a virtual 
photon carrying large three-momentum. Thus to obtain an adequate description, one must 
describe nucleons with relativistic wavefunctions, e.g., solutions to the Dirac equation. This 
in turn introduces the difficulty of treating antinucleon degrees of freedom, which we deal 
with approximately. Other issues also arise at higher energies, including the composite na- 
ture of nucleons themselves. This includes both nucleon electromagnetic form factors, which 
we include, and inelastic excitations of nucleons, which we do not. We also do not consider 
the effects of exchange currents among nucleons, e.g., as carried by charged mesons. These 



have been studied by Schiavilla et q?., |ll|| for example. For the purpose of this paper, the 
nucleus consists of interacting relativistic nucleons, considered to be elementary particles 
with form factors. 

A number of authors have considered relativistic extensions of the NRCSR with the 



approximations just mentioned. Walecka||12|| made the most direct extension by integrating 
the relativistic Coulomb response function over all energies to allow the use of closure. This 
results in a sum rule which is nearly identical in form to the NRCSR. MatsuijllJ showed, 
however, that such a sum rule would never be saturated by electron scattering, for which 
u>< |q|. He showed that, in the Fermi gas model, excitations of the Fermi sea (N scattering) 
are found entirely in the spacelike response, while excitations of the Dirac sea (NN pair- 
production) are found entirely in the timelike response. Unlike the NRCSR, for a system of 



Z Dirac protons the spacelike relativistic Coulomb sum tends to Z/2 in the high-|q| limit as 
a result of this separation. This work initiated interest in spacelike Coulomb sum rules based 
on energy transfers of uj < |q| only. Several authors have studied the case of a Fermi gas 
target. [13-16] We will see that the exclusion of timelike contributions from the Coulomb sum 
introduces kinematic factors related to relativistic nucleon recoil which must be accounted 
for in the formulation of sum rules for large |q|. 



DeForest|17| recognized the issue of relativistic nucleon recoil and proposed a correction 
in the form of a modified electric form factor: Ge(Q 2 ) = G e{Q 2 )\J (l + r ) / (1+2t), where 
r = Q 2 /4M 2 , which would be factored from the response function in the definition of the 
relativistic Coulomb sum. He conjectured that the resulting sum rule should have an inter- 
pretation similar to the NRCSR, but did not show this. DeForest's prescription has been 
used in the analysis of recent data.Jl8| However, we will show (in Sec. 5) that a result similar 
to his can be derived (RCSR-I) and does lead to a sum rule which bears close resemb- 
lence similar to the NRCSR, but is only accurate for Dirac nucleons, i.e., nucleons without 
anomalous magnetic moments. Donnelly et al.^E\ have also proposed a modification to the 



definition of the spacelike Coulomb sum. By requiring that their modified Coulomb sum 
tend to Z in the large- |q| limit, an integral equation for the modified form factor can be 
derived, then solved via an expansion in moments of the nucleon momentum. In lowest- 
order, this corresponds to the DeForest prescription evaluated at zero nucleon momentum. 
However, in higher order, this method introduces an energy dependence into the Coulomb 
sum function which then is not equivalent to a non- energy- weighted sum rule of the NRCSR 
form, from which correlation information can be extracted. We discuss these issues further 
in the concluding section. 

In this paper we derive a spacelike Coulomb sum rule valid for arbitrary three-momentum 
transfers, in which the relativity of the nucleons is taken into account. We assume that 
nucleons are the only degrees of freedom contributing to the spacelike sum and ignore, for 
example, antinucleon degrees of freedom, meson exchange currents and internal excitations of 
nucleons. In Section 2, we review the formalism for electron scattering from nuclear targets in 
the first-order Born approximation, including nucleon form factors with an assumption about 
their off-shell continuation. In Section 3, we consider elastic scattering from a single, free 
nucleon at rest as the simplest example of relativistic nucleon recoil effects in the constant-q 
response; this is shown to become a large effect for |q| ~M. In Section 4, we develop our 
"nucleons-only" approximation and derive a relativistic Coulomb sum rule which is exact 
in its treatment of nucleon recoil and Fermi motion. We give an explicit form for the two- 
body correlation function. In Section 5, we derive an expansion of the one-body term in 
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moments of the nucleoli momentum. This leads to practical sum rules for the analysis of 
experimental data. In Section 6, we present numerical results, including a test of our moment 
expansion in the Fermi gas model. Section 7 contains a discussion and conclusions, including 
a recommendation for the application of these sum rules to data analysis. 

2 Formalism for (e, e') Scattering 

We begin with the first-order Born approximation for the inelastic scattering of ultrarela- 
tivistic (|k| y> m e ) electrons from nuclear targets. Initial and final electron four-momenta 
are denoted by k^ 1 = (E, k) and k lil = (E', k'), and initial and final target four-momenta are 
denoted by P M = (P°,P) and P'^ l = ( y P' , P'), respectively. The differential cross-section can 
be written in factored form: 

d 2 a a 2 Ik' I 

' 'Lf(A;',fc)WV(P',P), (2.1) 



dn'dE' q 4 |k 

where g M = (u, q) = k^ — k nt = P" 1 — P M is the four-momentum transferred to the target in 
the laboratory frame by a single virtual photon, and a is the fine-structure constant. The 
lepton current tensor for an unpolarized electron beam is: 

L^(k',k) = X - £[Mk')7 M « s (k)] [uAk')Yu s (k)}\ (2.2) 

ss' 

where the fermion spinors u s (k) are given in (|A3|). 

All of the interesting target physics is contained in the nuclear response function: 

WV(P, q) = -£(f\J,(q)\i){f\X(q)\i)*8{u - E f + E t ), (2.3) 

/ 

where \i) and |/) denote initial^ and final nuclear many-body states, respectively. The 

tensor W^^P, q) may be decomposed into longitudinal and transverse response functions, 

with respect to the polarization of the virtual photon in a given frame: normally the target 

rest-frame. (We will not discuss the transverse response further in this paper.) If the z-axis 

is chosen along the three-momentum transfer q, i.e., q fl = (u),0,0, |q|), then the longitudinal 

polarization vector, chosen to obey e^-g = (Lorentz gauge) and e\ = 1, can be written 

e£ = (|q|, 0, 0,u)/y/Q^. The longitudinal response function is then defined 

W L (u;,q) = e1*W IMU e I i. (2.4) 



2 Here, as in Section 1, we assume a nondegenerate ground state \i). 



For theoretical reasons, it will be more convenient to work in the Coulomb gauge, defining 
the Coulomb response function 

Wc(u, q) = WooiuJ, q) = ^ W L (u, q), (2.5) 

where the second relation reflects this change of gauge. It is the Coulomb (not longitudinal) 
response function which appears in ( |1 . 1| ) . 

In general, the current operator J^{q) includes contributions from both charged nucleons 
and charged mesons in the target ground state, as discussed in the introduction. However, 
for the present work we ignore contributions to (|2.3| ) from meson currents. This allows us 
to write the nucleon current operator in the one-body form: 

J^q) = Jd 3 x e iqx ^(x)r^(x), (2.6) 

where ip(x) is the Dirac nucleon field operator in the Schrodinger picture. In general, the 
field operators and therefore the current operator ( |2.6| ) involve both nucleon and antinucleon 
degrees of freedom, as well as a sum over isospin projections. 

In the relativistic impulse approximation, the field operators in ( |2.6j ) refer to free particles. 
It is then natural to expand these field operators in terms of free Dirac wavef unctions, as in 
(|A2|). The 'yNN vertex operator T^ may then be specified in terms of its matrix elements 
between Dirac plane- wave spinors: ■u s /(p / )r / /u s (p), therefore it is conventional to express T^ 
in terms of Dirac matrices: 

r fi = F l ^ + i^F 2 a flu q", (2.7) 

where k is the nucleon anomalous magnetic moment and M is the nucleon mass. In many- 
body calculations involving both protons and neutrons, ( |2.7|) includes an implicit sum over 
proton and neutron isospin projections. In general, the form factors F\ and F 2 are scalar 
functions of p, p' and q. 

For scattering from free nucleons, the form factors i*\ and F 2 can be shown to be functions 
only of the scalar variable Q 2 = q 2 — uj 2 . These form factors are obtained from elastic eN 
scattering data, usually in terms of the more convenient electric and magnetic form factors: 



G E {Q 2 ) = F^Q 2 ) - ktF 2 (Q 2 ) 

G m (Q 2 ) = F^Q 2 ) + kF 2 (Q 2 ) (2.8) 

respectively, where t = Q 2 /4M 2 . Then F% and F 2 are given by 



2) = G E {Q 2 ) + rG M {Q 2 



f 2 (q 2 ) = G M {Q 2 )-G E m 

1 + T 

with the normalization convention that F[ (0) = i 7 f(0) = -F 2 "(0) = l and F"(0)=0. 

The continuation of these form factors for interacting nucleons is not unique, requiring 
dynamical information not contained in the impulse approximation. A common assumption 
is to use the on-shell values Fi(Q 2 ) and F 2 (Q 2 ) in constructing the vertex operator (|2.7|) . 



An equally reasonable assumption, which we adopt in this paper, is to evaluate only the 
electric and magnetic form factors at the off-shell variable Q 2 , i.e., the off-shell F\ and F2 
are defined: 



Fi(Q 2 ,Q 2 



G E (Q 2 ) + fG M (Q 2 ) 
1 + f 



F 2 (Q 2 ,Q 2 ) = G ^Q 2 )~ G ^Q 2 ) (2 . 10 ) 

1 +T 



where f = Q 2 /AM 2 and Q =q — (E p+Cl — E p ) 2 . Here E p = y/p 2 + M 2 is the energy of a free 
nucleon with three-momentum p. 

The crucial feature of the off-shell prescription fl2.10| ) is that the photon energy u enters 
r M in ( fT 7 ?!) only through the form factors G E (Q 2 ) and Gm{Q 2 )- Consider the limit of point 
nucleons, for which the form factors Fi(Q 2 ) and F 2 (Q 2 ) in ( |2.7| ) become constants: then the 
plane-wave matrix elements u s '(p+q)r^u s (p) depend only on the three-vectors p and q. For 
extended particles, both the F's and the G"s depend on the scalar Q 2 , and may also depend 
on another scalar, which we take to be Q 2 . (This choice is not unique, but is sensible given 
that p and q are the only other variables at the jNN vertex in the impulse approximation.) 
One may assume that either the F's or the G's, but not both sets, are functions of Q 2 only; 
we have made the second choice in fl2.10| ). This prescription is necessary to arrive at the 
particularly simple forms for the relativistic Coulomb sum rule which follows. 

We adopt the following parameterizations [fL9H for the nucleon electric and magnetic form 
factors appearing in ( p,10|) : 



(2.11) 
(2.12) 
(2.13) 
(2.14) 

where k p = +1.79 and K n —— 1.91 are the proton and neutron anomalous magnetic moments, 
respectively. Although the neutron electric form factor Ge,u{Q 2 ) is not identically zero in 
Ref. |fl9"| , it is known to be small compared to Ge, p (Q 2 )- Other nonzero forms for GE, n (Q 2 ) 



Ge, p (Q ) = 


= (l + Q 2 /.71GeV 2 )- 2 


Gm, p (Q ) = 


= (1 + k p )G e , p (Q 2 ) 


Gm,ti{Q ) - 


= K n GE,p\Q ) 


Ge,u{Q ) - 


= 



have also been given, fl9|, |20| and could be adopted for use in the present work. However, 
we shall make use of the fact that all form factors have Ge, p (Q 2 ) as a common factor to 
simplify the form of the relativistic Coulomb sum rule to follow. We note, however, that this 
proportionality cannot be exact for nonzero Ge,u{Q 2 )i since (jE in (0)=0. 

3 Single-nucleon Coulomb Scattering 

In this section, we consider the simplest nuclear target, one free nucleon, in the response 
function formalism of the previous section. This allows us to illustrate directly some of the 
features of relativity which also occur in the general nuclear case, and suggest how they may 
be treated in the formulation of relativistic Coulomb sum rules. In this case, the scattering is 
elastic on the target nucleon; the differential cross-section in the laboratory frame, which may 
be obtained from ( |2.1|) , is the Rosenbluth formula with form factors ( |2.8| ). We are interested 
only in the Coulomb scattering, which is expressed in terms of the Coulomb response function 
Wc(uJ,q) = W 00 (uJ,q). 

For a free nucleon, the Coulomb response function ( |2.5|) may be evaluated from ( |2.3|) and 
( p.6|) by specifying that initial and final nuclear states are simply Dirac plane-wave states of 
momenta p and p' = p+q, respectively. Using the nucleon field operator of (|A~^) , we find: 

Wg(p, q) = J^ S(u-E p+Cl + E p ), (3.1) 

where the 00-component of the nucleon current tensor is given by 

£oo(P>?) = 2 H[^'(p') r o« s (p)][Ms'(p / ) r oM s (p)] , (3.2) 

ss' 
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and the the energy factors in the denominator of ( |3.1| ) reflect our convention that the fermion 
spinors ( |A3| ) are normalized to 2E particles per unit volume. The delta-function in ( |3 . 1| ) 
guarantees energy conservation for the free nucleon states: u = E p+(i — E p . 

Notice that only nucleon, i.e., not antinucleon, states appear in ( |3.2| ). The energy required 
to excite a free particle of momentum p into a state with momentum p+q satisfies -E p+q — E p < 
|q|, hence the process N ^N is restricted entirely to the spacelike region uo< |q|. In contrast, 
the energy required to excite an NN pair from the vacuum satisfies uj = E p+(l +E p > |q|, and 
hence is restricted entirely to the timelike region u > |q|. (Pair-production arising due to 
scattering from a free nucleon may be considered either an excitation of the nucleon itself 
or a two-body, i.e., meson exchange, process and is ignored here.) Thus the antinucleon 
components of the field operators of ( |A2|) do not appear in the (spacelike) response function 



for a free nucleon. We shall find it advantageous to use this separation of excitation spectra 
to explicitly remove the contribution of NN pairs from the relativistic sum rule. (This issue 
was first raised by Matsui.[|T3j]) 

Returning to the Coulomb response function, we first give an explicit form for the 00- 
component of the nucleon current tensor (|3.2|). After substituting ( |2.7| ) and ( [2.10| ) into (|3.2| ), 
a straightforward trace calculation yields 



L&P, Q) = ^^(E p+q +E p f + ^^1 [f (£ P+q +£ P ) 2 - (1 + f )q 2 ] , (3.3) 

1 + r 1 + r L J 

where, of course, for a free nucleon f — r and Q 2 = Q 2 . If we put the initial nucleon at rest, 



evaluation of (|3.3|) at p = yields 



£&(P.?)l P =o= 2M ^ + M ) G e{Q*), (3-4) 

where we have used f = (E^ — M)/2M at p = 0. Thus stationary nucleons couple only 
through their electric form factor. (A sensible feature of the off-shell prescription (|2.10| ) is 
that this reduction also occurs for interacting nucleons.) The response function is obtained 
by substituting ([0|) into flnp. 

The Coulomb sum function is obtained by integrating Wc(u>, q)/G 2 E JQ 2 ) over energy 
transfer u, as in ( |1.5| ). For electron scattering from a free nucleon, it is natural to restrict 
the integral to include only spacelike (u><|q|) photons, thus excluding NN pair production: 



il , Wgfaq) 



duo 



GUQ 2 



E^ + M __ M q : 
p =o 2_E q Eq Q 



(3.5) 
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where in the second relation we have made use of the nucleon energy-momentum relation 
at p = to express the result in a form which lends itself easily to interpretation. The first 



factor M/Eq has its origin in the energy denominator -E p + q in ( |3.1| ), which in turn comes 
from the relativistic normalization of the fermion spinors (|A3|). The second factor qf/Q 2 
arises from our choice to work in Coulomb gauge, as can be seen by comparing with ( |2.5|) . 
Taken together, these factors represent the effect of relativistic recoil (in Coulomb gauge) 
for a single nucleon. The first relation in (|3.5|) shows that this recoil effect tends to unity 
in the nonrelativistic limit |q| <kM, and tends to the value 1/2 in the ultrarelativistic limit 
|q| y> M. (It is in this sense that Coulomb gauge is a convenient choice, since the analog of 
( |3.5|) in Lorentz gauge tends to zero in the limit |q| y> M.) 

Matsui|l3| and DoDang et a/.|l4| have previously noticed this limiting value of 1/2 for 



the spacelike Coulomb sum, for a Fermi gas of Dirac (Ge = Gm — 1) nucleons. The Fermi gas 
result is easily obtained by inserting ( |3.3| ) into ( p.l| ) and integrating over the Fermi sphere 
|p| <Pf with an appropriate factor to account for Pauli blocking, as is done in ( p.2| ). It is 
easy to show that this quantity approaches Z/2 (for Z protons) as |q| — ^oo. That this is a 
purely kinematic effect of relativistic nucleon recoil is clear from the single-nucleon example 

The high-|q| behavior of ( |3.5| ) is in contrast to that of the NRCSR ( |1 . 1 1| ) , for which 



Snr((1) —> Z (Z = \ here) as |q| — »oo, suggesting that in the relativistic limit, the Coulomb 
sum does not directly count the number of charged scatterers, as in the nonrelativistic case. 
A simple modification of ( |3.5| ) would, however, seem to remedy this for the case of a single 
nucleon: divide both sides by the kinematic recoil factor on the right-hand side, so that 
the large- 1 q| limit is now unity, as would be expected for a sum rule. That this actually 
corresponds to a sum rule for the many-body case, with a suitable modification for the 
nucleon momentum distribution (Fermi motion) in the target, is shown in the next two 
sections. 

A final note: For a free nucleon, the delta-function of ( |3.1|) appears in all components 
of W^ u . To obtain the Rosenbluth formula, one must integrate fl2.1|) over dE' (or du) at 
fixed electron scattering angle 8, rather than at fixed q, as in the Coulomb sum (|3.5| ). The 
constraint that to and q are not independent variables, but are related by Q 2 = 4EE' sin 2 (9/ 2), 
then gives rise to the familiar recoil factor: 



E~ 



2E . „/0 



1 H sin 



(3.6) 



M V2, 
which reduces to unity in the nonrelativistic limit E « M, where E is the electron beam 

energy. Thus it is well-known that not only form factors, but also nucleon recoil factors, 
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modify the differential cross-section for electron scattering in the relativistic domain. 

4 Relativistic Coulomb Sum Rule for Interacting Nu- 
cleons 

In this section, we derive a sum rule for a target of interacting nucleons, where the relativity of 
the nucleons is taken into account. Following the example of a single nucleon in the previous 
section, we find that if we restrict the sum to spacelike four-momentum transfers, which 
corresponds to (e, e') scattering experiments, we are led to a "nucleons-only" approximation 
in which antinucleon degrees of freedom are ignored. The new feature in the sum rule given 
here is the kinematic correction for relativistic recoil, which depends on the momentum of 
the nucleons in the target. In the following section, we show how this correction can be 
calculated in successive approximations in terms of moments of the nucleon momentum. 

We begin with the Coulomb response function in the rest-frame of a many-body nuclear 
system, which is given by W 00 (P,q) of ( |2.3| ), evaluated at P = 0: 



Wc{w, q) = £ \M fi (u>, q)| 2 S(u -E f + E t ), (4.1) 

/ 



where the Coulomb transition matrix element is: 



M fi (u,q) = (/|Jo(w,q)|*> = J d 6 x e^ x (/|V(x)r ^(x)|z). (4.2) 

These are the relativistic equivalents of ( |1.3|) and ( |1.4j) . The jNN vertex operator T is 
given by ( [2.71) an d includes an implicit sum over isospin projections, and the field operator 
V>(x) is given in ( |A2^ ) and includes both nucleon and antinucleon degrees of freedom. As 
discussed in Section 2, the nucleon form factors are assumed to be given by ( |2 . 1 0|) such 



that the only u dependence in Jo(q) enters through the nucleon electric and magnetic form 
factors, Ge{Q 2 ) and Gm(Q 2 )- However, given the specific functional forms in ( |2.11| )-( p.l4j ), 
one can instead consider all the ui dependence in ( |4.2j ) to enter only through the common 
form factor Ge, p {Q 2 )- Dividing ( |4.2j ) by Ge, p (Q 2 ) therefore removes all the uj dependence 
from the nuclear matrix element M/j. Dividing ( |4.1| ) by G E „{Q 2 ) then allows us to generate 
a non-energy-weighted sum rule, as follows. 

We consider the Coulomb sum function defined by 

ylql W c {u>, g) ,. _. 
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where the integration is performed only over energies u < |q|, corresponding to electron 
scattering experiments, as in ( |3.5| ). This is unlike the nonrelativistic Coulomb sum ( |1.5| ), in 
which the integration is over all u> > us e i in order to invoke closure. Substituting ( }4.1| ) into 
(|4.3|) and performing the integration, we obtain: 



jf:w<|q| 

s(q)= E 

ffr 



M fi {u,q) 



2 



GeAQ 2 ) ' (4 ' 4) 

where the sum over all final states in (|4.1|) is now reduced to a sum over only those final 
states / 7^ i which are accessed by spacelike four-momentum transfers. For the case of a free 



nucleon target at rest, (|4.4j) reduces to (|3.5|) . 

The nucleon field operator ( |A2|) , and therefore the transition matrix elements (|4.2| ), in- 
volve both nucleon and antinucleon degrees of freedom. However, because electron scattering 
experiments produce spacelike virtual photons exclusively, the resulting nuclear excitations 
are restricted to those whose excitation energies satisfy u = Ef — Ej < |q|, as indicated in 



We have already seen in (|3.5|), for the case of a free nucleon, that restricting the energy 



transfer oj removes the contribution of NN pair production from the Coulomb sum, which 
would enter if u>> |q| were included in ( f4.3p . Similarly, for a noninteracting Fermi gas, only 
excitations involving a p , s ,a ps (N scattering) enter into the spacelike sum. This is because 
there are no antinucleons present in the renormalized nuclear ground state, thus the terms 
involving b ps b p / s / (N scattering) and b p ' s 'a ps (NN pair-annihilation) vanish identically, and 
because the terms involving a p , s ,b ps (NN pair-creation) are restricted entirely to the timelike 
region, since u = E p+q +E p > |q|. Furthermore, the separation of nucleon and antinucleon de- 
grees of freedom is complete in the Fermi gas model: not only are antinucleon scattering and 
NN pair terms completely removed from the spacelike sum, but the N — > N contribution is 
guaranteed to lie completely in the spacelike region through the relation u = E p+Ci —E p < |q|. 
(This feature is unaffected by the inclusion of scalar and vector potentials like those in Quan- 
tum Hadrodynamics, for example, as noted by Do Dang et a/.0]) Thus, in the Fermi gas 
model, the spacelike nuclear response (}4.4|) can be represented by effectively dropping all 
reference to antinucleon degrees of freedom at the level of the field operator (|A2|) . 

For an interacting nuclear system, however, the separation of nucleon and antinucleon 
degrees of freedom is no longer exact. In general, both N and N degrees of freedom can 
contribute to the spacelike nuclear response, and some of the nucleons-only response N ^N, 
which is entirely spacelike in the Fermi gas model, can be pushed into the timelike region by 
two-particle interactions. However, we can assume that true NN pair production will still 
occur mostly for lo> |q|, and therefore can be neglected in the Coulomb sum ( [4.4| ). Although 
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tests of this assumption are outside the scope of this paper, we note that it is satisfied 
exactly for infinite nuclear matter in the RPA, i.e., a spacelike photon excites only spacelike 
excitations even through the summation of an infinite number of RPA ring diagrams, and 
does not produce NN pairs. Antinucleon contributions, i.e., those involving frL&pv and 
b p i s ia ps , enter as relativistic effects in the nuclear target structure. In Hartree calculations 
for finite nuclei, these terms make only small (~ 1%) contributions to the relevant (vector) 



density, [^TJ and shall also be neglected in this paper. Therefore, we adopt a "nucleons-only" 
approximation, in which antinucleon degrees of freedom are removed from discussion at the 
level of the field operator ( |A~2f) . In effect, we are assuming that the dominant relativistic 
effect comes from single-nucleon recoil at high |q|, and not from relativistic aspects of nuclear 
structure. Thus the initial and final nuclear states in (|4.4j ) will be treated as interacting 
many-nucleon states without N components. 

Substituting the field operators given by (|A2|) into the transition matrix element 



we have in our "nucleons-only" approximation: 

M /i (w,q) ^G E:P (Q 2 )J2J2js'sAP^)(f\ a i+ n s'a a Psa\i), (4.5) 



p+qs' 



where the "reduced" spinor matrix element is defined: 



, _ ?V CT (p + q) r (p,g) UsAjp) , v 

Pq ' l2E^ G E>p mj2E p [ } 



In (|4.5| ) and ( [4.6|) , nucleon spin states are represented by Latin indices and nucleon isospin 
states are represented by Greek indices. The operator r (p, q) is diagonal in isospin space, 
therefore the transition matrix element ( |4.5|) simply involves a sum over proton and neutron 
isospin projections. The term "reduced" in (f4.6|) refers to the division by Ge, p {Q 2 ), following 
( |4.3|) . The matrix element ft4.6|) is then a function only of p and q, i.e., not of u, through 



our assumption of proportional form factors. Substituting (|4.5|) into (f4.4|) leads to: 



E ( f i) = Y.J2J2J2 Jr'r, P (p' 5 q) J S 's,<x(P,q) 

pp' rr' ss' &P 

{ J2 ( i l a P 'rpV+qr'pl/)(/l a p+q S 'a a p^K) 

— (i\a plr a p i +cir >p\i)(i\a p+/cr ap S<7 \i) >, (4.7) 
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where we have subtracted the elastic term f = i explicitly, in order to obtain a sum over the 
complete set of nuclear final states. For notational simplicity, we assume a nondegenerate 



ground state \i). Equation (477) is then the most general form for the Coulomb sum function 
( |4.3| ), given our nucleons-only approximation and our assumptions about nucleon form factors 
and their off-shell continuation. 

In order to obtain a sum rule from ( |4.7|) , we must sum over a complete set of final states 
\f) and use closure. The completeness of such a set necessarily includes excited nuclear states 
accessed by both spacelike (uj < |q|) and timelike (u> |q|) photons, which can be expressed 



E I/X/I+ E l/X/l 

/:uK|q| /:^>|q| 



(U 



It is consistent with the nucleons-only approximation imposed in (|4.5| ) to neglect the con- 
tribution of timelike states / :u> |q| to the closure of the sum in ( |4.7|) , since these states 
predominantly involve antinucleons. For a uniform Fermi gas in the Hartree-Fock approx- 
imation, we have already seen that the separation between N and N degrees of freedom 
guarantees that the spacelike Coulomb sum be completely exhausted in the nucleons-only 
approximation. Similarly, for an interacting Fermi gas in the RPA, nucleons are excited only 
for oo < |q|, and NN pairs are produced only for u> > |q|, thus the same argument holds. It is 
possible, however, that restricting the sum in ( [4.4] ) to spacelike states /:u><|q| does not in 
general exhaust the nucleons-only excitation spectrum, e.g., finite size effects and two-body 
interactions beyond the RPA may push some of the nucleons-only response into the timelike 
region. We are aware of no general argument which applies. We will simply assume that the 
nucleons-only excitation spectrum is saturated by states / : u> < |q|, so that we can invoke 
closure in ([4.7|) . Extension of the sum to include nucleonic excitations with u > |q| could be 



done by theoretical means, as in Schiavilla et a/.,|22| for example. Applying closure to the 
first term in ( f4.7|) removes the sum over final states, and the expression in curly brackets 
becomes: 



— \*l a p'rp a p'+q»*'pr/ \^l a p+qs'o- a pso-| V 



(4.9) 



Equation ( [4.7|) is now a sum rule, equating the sum (|4.3| ) to the ground state expectation 
values in 
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To cast the sum rule in a more transparent form, we separate out the one- and two-body 



parts of the operator in the first term of (|4.9|), using the momentum-space anticommutation 



relations for the nucleon creation and destruction operators: 

{apsa, flpVa'} = Spp'tiss'ScTCT', (4-10) 

where the anticommutator of any two creation or destruction operators vanishes. Moving 
all creation operators to the left in the first term of (f4.9|), we obtain: 



(i\a psa a psa \i)5 pp >5 r > s >5 pa 5 rs 

+ (*l a p+qs'o- a p'rp a p'+qr-'p a ps°-P/ 



V'\ a p'rp a p'+qr'p 



iKilop+WWl*)}' ( 4 - n ) 



where we have extracted S rs from the first matrix element, which is diagonal in spin projec- 
tion. Inserting ( [4.1 1|) into (|4.7|) allows us to separate the sum rule expression into one-body, 
two-body and elastic amplitudes, in complete analogy to the nonrelativistic expression ( |1.8|) . 
Here, however, we have used momentum operators, rather than local field operators, in order 
to eliminate the antinucleon components. 
We separate the sum rule as follows: 

E(q) = £«(q) + E( 2 )(q). (4.12) 

The first term gives the one-body contribution: 

£«(q) = 2]>> CT (p)r CT (p,q), (4.13) 

per 

where we have defined the momentum distribution function n ff (p) = (i\a psa a psa \i) for isospin 
projection a (which we assume is independent of s), and the relativistic recoil factor 



?-v(p,q) = ^X^'^P' 1 !) 



I L° (p,q) (4M) 



17 



where Lq (p, g) is the nucleoli tensor component ( |3.3| ). In the nonrelativistic limit |q|<CM, 
r CT (p, q) — » 1 and E^(q) —* Z, as in ( P-.8| ). For relativistic momenta |q| ^> M, however, we 
have rv(p, q) < 1, which represents the kinematic effect of relativistic nucleon recoil, as in 
( p.5|) . If the Fermi momentum in ( |4.13|) could be neglected entirely, i.e., if we were to set 
p = in ( f4.14|) , then EW(q)/Z would be given by ( |3.5| ). The second term of (|4.12|) is given 
by 



s(2) (q) = EEEEjrv>'q)j s v(p.q) 

pp' ss' rr' °"P 

\ V 1 1 a p+qs'cr a p'rp a p'+q>-'p a pso- r ) 



~~ (* I a p'rp a p'+qr'p I V V>\ a p+qs'o- a pso- 1 V j • 



This expression contains the full two-body probability in momentum space, in analogy with 
the second term of (|1.8|), but is weighted by the kinematic factor j*i rtP (p', q) j s 's,a(p, q)- 

To separate out the correlation part of E^ 2 ^(q) as in (|1.9[), we first evaluate ( |4.15| ) for an 
uncorrelated target ground state, for which we obtain 



s in(q) = -£tt|£X! j^,-(p,q)(*l a p+ qs v a P^K) 



TV ■ 

o- 1V °" p ss 



(4.16) 



where N„ = Z,N for protons and neutrons, respectively. (In the uncorrelated case, there is 
no contribution in (|4.16|) from cr^p.) We may then define the correlation function 



C(q)^E( 2 )(q)-E^)(q) 



(4.17) 



to give a measure of true two-body correlations, and rewrite ( 4.12 ) 



E(q) = E«(q) + C(q) + E( 2 )(q), (4.18) 

in analogy with ( |1 . 1 1| ) . This result will be referred to as the relativistic Coulomb sum rule 
(RCSR). 

The uncorrelated term E^ (q) is related to the square of the elastic target form factor, 
with the difference that the proton and neutron contributions are added incoherently in 
( |4.16| ). If, however, the contribution of neutrons to ( 4.16 ) can be neglected, then 



Sg(q) 



1 |^(q)| 2 

% GE,p{Qel) 



(4.19) 



completing the analogy of ( |4.18| ) to ( |1 . 1 1| ) . For example, if the contribution of the magnetic 
{Gm{Q 2 )) terms is zero, as discussed by FriarpBI for a spin-saturated target, then with 



Ge,u = 0, as assumed in ( |2.14|) , ( |4.19|) is exact. For Ge,„7^0, ( |4.19|) can be suitably modified 



under reasonable assumptions, e.g., similar neutron and proton distributions in ( }4.16|) . 

The RCSR given here is exact in its treatment of nucleon Fermi motion and in principle 
allows the identification of nucleon- nucleon correlations to arbitrarily high |q|, given our 
assumptions about antinucleon and meson degrees of freedom, nucleon excitations, and off- 
shell continuation of nucleon form factors. From ( |4.18|) , the extraction of the correlation 
function C(q) from the calculated Coulomb sum (|4.3| ) requires the removal of the nuclear 
elastic form factor, as just discussed, as well as a reliable evaluation of the one-body term 
( |4.13| ). For the latter, this form of the RCSR presupposes that one has in hand the nucleon 
momentum distribution n CT (p). 

5 Expansion in Moments of Nucleon Momentum 

The momentum distribution n CT (p) is not generally available directly from experimental data. 
In order to use the RCSR as given in the previous section, therefore, one would have to rely 
on a theoretical model to calculate the one-body term E(q). An alternative approach, 
which we develop in this section, is to expand the one-body term ( |4.13|) in moments of 
the nucleon momentum, i.e., in averages of powers of the nucleon momentum, weighted by 
the distribution n CT (p). The one-body term ( |4.13|) is then replaced by a sum in terms of 



these moments, the first few of which may be known. This procedure leads to a series of 
approximate sum rules, each of which depends on higher momentum moments. In each case, 
it is possible to identify appropriate modifications to the definition of the relativistic Coulomb 
sum (|4.3j ) and arrive at a series of modified sum rules, which share certain features with 
the NRCSR. Of course, such a procedure is useful only if the moment expansion converges 
rapidly, and if an adequate number of moments can be obtained reliably. In the next section, 
we will demonstrate that this expansion indeed converges quickly for a uniform Fermi gas, 
and argue in the following discussion that this is likely to be the case in general. 

5.1 Sum Rule I 

The lowest-order momentum moment expansion is obtained by setting p = in ( |4.14|) : 



r ff (p, q) ~ tv(0, q) = -^= S ap , (5.1) 



2£ q 
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which is just the factor that appeared in ( |3.5| ) for a single proton at rest. This is equivalent 
to approximating the reduced matrix element (f4.6|) by its value at p = 0: 



is' S)CT (p,q) ^ js>sA°> q) 



l Eg + M 

ss i0 ap . 



(5.2) 



Since we evaluate at p = 0, ( pM| ) and (|5.2[ ) do not involve magnetic terms. The absence of 
neutron electric contributions is a result of our approximation Q2.14 ). 
Inserting ( |5.1| ) into ( |4 . 1 3| ) , the one-body term becomes 



E«i 



q) 



E a + M 



2E n 



(5.3) 



Similarly, inserting ( |5.2| ) into ( |4.15| ) and (|4.16| ), the two-body correlation term ( |4.17|) becomes 



C(q) 



E n + M 
2E n 



2-^2.^} (i\ a p+(is a p > s > a p'+cis'aps\i) 

pp' ss' 



\*l 0, p's' 0, p'+qs'r/\' i l c, 'p+qs a 'psr/ 



(5.4) 



where the proton isospin label has been suppressed for clarity. 

Since both E^q) and C(q) now appear with the same overall kinematic factor, i.e., that 
which appeared in (|3.5| ), which is a function only of q, it is possible to modify the definition 
of the Coulomb sum (fOf) and obtain a more conventional form for the sum rule. We define 
the modified Coulomb sum function: 



SV(q) = 



E(q) 



mqj r^q) 
where the lowest-order recoil correction factor is 



\l\ W c (uJ, q) 

■5 GIAQ 2 ) 



(5.5) 



n{q) 



E n + M 
2E n ' 



(5.6) 



and Wc(&, q) is to be extracted from experimental data, using (|1.1|). We then find the 
approximate sum rule: 



&(q)~Z + CIr(q) 



1 



r/(q) 



i l^(q) 



^ GeAQ 



el) 



(5.7) 



This result would be exact for a target of free, stationary protons. The two-proton correlation 
function in ( |5.7| ) is given by 
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Cl(q) 



pp' ss' 



''p+qs^p's' 



7 j 7 , 1 \2| a p+qs a p's' a p'+qs' a p.sl2 ; 



Z-l 



^ a. 



p'.s 



/Op' + 



qs '|i)(i|4 



p+qs a p.' 
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in close analogy with (|1.9| ). Expression ( |5.7| ) will be referred to as RCSR-I. 

Like the NRCSR given in ( |1 . 1 1|) , RCSR-I involves a one-body term Z which simply counts 
the total number of charged scatterers, and a two-body term which contains information on 
nucleon-nucleon correlations. The function C/(q), however, is not the Fourier transform of 
the usual spatial correlation function, as defined in terms of local nucleon field operators, 



i.e., simply inserting the Dirac nucleon field operators of (|A2|) into (|1.9| ) will not yield (|5.8|) , 
because of the antinucleon degrees of freedom which are eliminated from ( |5.8|) , and the 
kinematic factors related to relativistic nucleon recoil which result from that elimination. 
However, if the nonrelativistic result ( |1.9|) is expressed in momentum space, it is exactly of 
the form ( |5.8|) . Thus this version of the RCSR, although based on a rather strong assumption 
about Fermi motion in the target, is closest in form to the familiar NRCSR. 



5.2 Sum Rule II 

The procedure which led to RCSR-I can be extended to include higher moments of the 
nucleon momentum. In this subsection, we keep terms only through second order. We first 
expand the relativistic recoil factor: 



rv(p, q) = rl + r^ (p ■ q) + r 2 CT p 2 + r 3 CT (p ■ q) 2 + (9(p 3 ), (5.9) 

for isospin projection a. Since the recoil factor ( |4.14j) is a function only of p and q, i.e., not 
u, the expansion coefficients r" are functions only of q. Explicit forms for the r" are given 
in Appendix B, where we find rj = r™ = as a consequence of approximation ( |2.14|) . 
Substituting (|5\9| ) into the one-body term (|4.13|) leads to 



Eg>(q) ~ 2£ n p (p)[rg+ (rg + ^rf)p 2 ] + n n { V )[(r n 2 + \r^ 



q- 

3 



(5.10) 



where we have used the spherical symmetry of n a (p) to eliminate the r^ term and integrate 
the angles: (p • q) 2 — >|p 2 q 2 - We car 
i.e., the second momentum moment: 



the angles: (p • q) 2 — >|p 2 q 2 - We can now identify the average squared nucleon momentum 
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(p\ = 



N n 



2]>> (T (p)p 2 



(5.11) 



where N a = Z, N for protons and neutrons, respectively. The one-body term ( |5.10| ) can then 
be expressed 



-(i) 



E^(q) 



I + {A + £■§) (p 2 ) J +N\(r n 2 + £•?) ( P 2 )„ 



(5.12) 



We thus have an approximate expression for the one-body term ( |4.13| ) which depends only 
on the two parameters (p 2 )cr, and reduces to 



if(p 2 ) P =(p 2 )n = 0. 



We can now follow a procedure similar to that of ( |5.5| ), by defining the modified Coulomb 
sum function: 



S 



//(qj 



S(q) 



M du Wb(w,q) 



where the second-order relativistic recoil factor is 



Gi„(Qiy 



r E,p\ 



(5.13) 



r/z(q) = 



,2 N . „2 



-§ + (ri + V") (p 2 )p + y (rj + \ri) (p\ 



We thus find the second-order approximate sum rule: 



S'ii(q)^ + C'j/(q)- 



rj/(qj 



i l^(q)l 2 

^ ^E,p(QelJ 



(5.14) 



(5.15) 



in analogy with 



The second-order nucleon-nucleon correlation function is given by 



C7/(q) 



c(q) 



(5.16) 



where C(q) is defined in ( |4.17| ). Expression (|5.15|) will be referred to as RCSR-II. Explicit 
expressions for the two-body term S(q) and the correlation function C(q) could be given to 
the same order of approximation made in ( |5.9| ), but would require the expansion of j s > s ^(p, q) 
to second-order in p. Since we shall not use these explicit forms, we do not pursue this here. 
Although RCSR-II enjoys a less direct correspondence with the NRCSR than does RCSR- 
I, it is more accurate in its treatment of Fermi motion. More importantly, as we will see 
in the next section, the higher-order terms which are included in RCSR-II introduce effects 
from proton and neutron anomalous magnetic moments which are appreciable. 
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6 Numerical Results and Test in Fermi Gas Model 

In this Section, we present numerical results for our relativistic Coulomb sum rules. We first 
consider the recoil correction factors r/(q) and rjj(q), which appear in RCSR-I and RCSR-II, 
respectively, and find that terms ~ (p 2 )<r in our momentum expansion bring in anomalous 
magnetic moment effects which are appreciable. We then illustrate in the Fermi gas model 
the accuracy with which RCSR-I and RCSR-II allow the evaluation of the one-body term. 

Ideally, one would like to extract (p 2 ) CT directly from experimental data. Since that is not 
always possible, however, one may have to rely upon estimates based on theoretical models. 
The simplest such model is a uniform Fermi gas, where n a (p) = 9(p F — |p|). In that case, we 
have: 

(P 2 ). = | Pf\ (6-1) 

where p F is the Fermi momentum for isospin projection a. In this example, we assume Z = N 
and take p F = pp = 1.42 fm _1 = .28 GeV/c to match the average density in the interior of 
finite nuclei. This gives (p 2 ) a — -047 GeV 2 /c 2 . For Z ^ N, it is reasonable to assume that 
the Fermi momenta for protons and neutrons scale according to p f =pf(2N ct /A) 1 ' 3 . 

Figure 1 shows the recoil correction factors r/(q) and rjj(q), given in ( |5.6| ) and ( 5.14 ), as 



functions of the three- momentum transfer q. The dotted curve is 77 (q), and is indistinguish- 
able (on the scale of this plot) from the factor given by DeForest: G 2 Ep /G Ep = (1+t)/(1+2t), 
suitably averaged over angle. The dot-dashed curve is rv/(q) with k p = k u = 0, i.e., for Dirac 
nucleons. Both curves tend to the value 1/2 as |q| — * oo. Their difference is a measure of 
the importance of second-order terms in the moment expansion when anomalous magnetic 
moments are ignored, and is at most ~l/2% anywhere in the range 0< |q| <oo. The dashed 
curve is rjj(q) with k p = 1.79 and n n = 0, and shows that the proton anomalous magnetic 
moment leads to a noticeable enhancement in the recoil factor: 5% at |q| = 1 GeV, and 
reaching 24% as |q| ^oo, when compared to rj(q). The solid curve is r//(q) with k p = 1.79 
and n n = — 1.91, and shows a similar enhancement due to the neutron anomalous magnetic 
moment when N = Z: 2% at |q| = 1 GeV, and reaching 13% as |q| — >oo, when compared to 
r 7 (q). 

It is clear from Figure 1 that rjj(q) is significantly different from 77 (q) only when anoma- 
lous magnetic moments are included. This can be understood by considering the coefficient 
of (p 2 ) p in ( |5.14| ), i.e., the quantity rf+ (q 2 /3)rf, where rf and rf are given by flB3"|) and 



([ETD, respectively. We find that the coefficients of (1 + k p )° and (1 + k p ) 2 in that quantity 
are of roughly the same magnitude and have opposite sign. Thus the second-order pro- 
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ton contribution to ( |5.14 ) is nearly zero when k p = 0, and rjj(q) ~ T/(q). However, when 
K p — 1.79, we have (l + /t p ) 2 ~8. In that case, the coefficients of (1 + /%,)° and (1 + k p ) 2 no 
longer cancel, and rjj(q) is significantly enhanced. Similarly, for neutrons, both r^ and r% 
are positive and proportional to /t 2 . Thus when K n = —1.91, we have k 2 ~ 4 and r/j(q) is 
further enhanced. From the specific functional form of Lq (j>, q), as given in ( |3.4| ), we note 
that terms in ( |5.9|) which are of higher order in p 2 will not involve higher powers of the 
anomalous terms (1 + k p ) 2 and k 2 . Therefore, we expect that the moment expansion will 
converge if the moments themselves converge. 

We now present a numerical test of our moment expansion in the simple system which 
includes Fermi motion: a uniform Fermi gas. The spacelike Coulomb response function for 
this system has been studied by several authors. [13-16, 24] Substituting (|3.3| ) into (|3.1|) , and 
integrating over the Fermi sphere |p| <Pf with an appropriate factor to account for Pauli 
blocking, we obtain: 

W c (u, q ) = 2 ^/(^^(ri^^tP+qll^^-^ + ^p)' (6-2) 



where the occupancy function n a (p) = 6(pp — |p|). Inserting ( |6.2|) into ( [4. 3D gives the rela- 
tivistic Coulomb sum: 

s (q) = 2 E/ 77m^(p)[ 1_ ^(p+ c i)] r ^(P' c i)' ( 6 - 3 ) 



where we have used definition ( |4.14j ). Comparison with Q4.13 ) makes the identification of 



one- and two-body contributions trivial for this model. Furthermore, the two-body term 
involves only Pauli correlations, which vanish for |q| > 2pp, and F e i(q) = for |q| > 0. We 
now demonstrate in several cases the accuracy to which S(q) of ( |6.3D is represented by 
RCSR-I and RCSR-II. We take (|6.3| ) to represent the experimental Coulomb sum obtained 



from (|4.3j ) and ( |1.1|) , and divide out the recoil factors r/(q) and r n (q), as in (|5.5|) and (|5.13| ). 



Figure 2 shows several versions of the Coulomb sum (divided by Z) for a free Dirac 
proton gas, i.e., k p = n n = 0, versus three- momentum q in units of the Fermi momentum 
Pf- The dot-dashed curve is the unmodified Coulomb sum E(q)/Z, given by (|6.3|) ; this is 
exactly the result given by Matsui,[]T3[ and tends to the value 1/2 as |q| ^oo. The dashed 
curve is Sj(q)/Z, as defined in (|5.5|), and the solid curve is 5 , //(q)/Z, as defined in ( |5.13| ). 



Here, as in Figure 1, the result given by DeForest is indistinguishable on the scale of this 
plot from RCSR-I. Both the dashed (RCSR-I) and solid (RCSR-II) curves approach unity 
quickly as |q| — > oo, although RCSR-I overshoots slightly (~ 1/10%) in the 1 GeV region. The 
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convergence to unity as |q| — >oo is guaranteed in the Dirac case, since both the unmodified 
Coulomb sum E(q)/Z and the recoil factors r/(q) and r/j(q) approach the value 1/2, for 

Kp = K n = U. 

Figure 3 shows the same quantities for a free proton gas, with k p = 1.79. In this case, 
E(q)/Z — » .62 as |q| — * oo, due to the proton anomalous magnetic moment. In contrast to 
the Dirac case, RCSR-I fails to level off at |q| =2pp, but rather Si(q)/Z — >1.23 as |q| — >oo. 
However, since r/j(q) accounts for the anomalous moment, RCSR-II levels off immediately 
at |q| = 2pp and Sn(c\)/Z — > .99 as |q| — > oo. The difference from unity is a measure of the 
importance of higher-order terms, e.g., (p 4 )<r, which have been neglected in our expansion. 
Figure 4 shows the same results for a symmetric (N = Z) Fermi gas. In this case, neutron 
terms further enhance the unmodified sum: S(q)/Z^.68 and Si(q)/Z — > 1.35 as |q| — >oo. 
However, when the neutron anomalous moment is incorporated into the recoil factor rjj(q), 
we again have Sjj(q)/Z ^ .99 as |q|— >oo. 

We have shown that it is possible to estimate reliably the one-body term in the Coulomb 
sum if one accounts correctly for nucleon recoil. Due to the simplicity of the Fermi gas 
model, however, the numerical study given here does not test the validity of our nucleons- 
only approximation, since the separation of iV and N degrees of freedom is exact in this case, 
nor does it test our assumptions about the off-shell continuation of nucleon form factors, 
since r = f in this case. This study does, however, give an indication of the importance of 
anomalous magnetic moment effects, and how accurate these approximate sum rules can be 
for the extraction of NN correlation information from electron scattering data, when the 
momentum distribution of the nuclear target is like that of a Fermi gas. 

7 Discussion and Conclusions 

The main result of this paper is the derivation of a Coulomb sum rule applicable to the anal- 
ysis of inelastic electron scattering experiments on nuclei at high three-momentum transfers. 
The restriction to spacelike four- momenta (u< |q|), appropriate for electron scattering, ap- 
proximately eliminates A^ A" pair-production from the response function. We then ignore the 
small contributions which arise due to antinucleons in the target ground state; this is our 
"nucleons-only" approximation. The sum rule is given in three forms. 

The most general form of the relativistic Coulomb sum rule (RCSR) appears in ( 4.18| ). 



As for the well-known non- relativistic case (NRCSR), the RCSR can be decomposed into a 
one-body term (£W(q)) and a two-body correlation term (C(q)) after removing the elastic 
term. The correlation term is of considerable interest, but is expected to be dominated by 
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the one-body term for all but the smallest three-momentum transfers. In the NRCSR, the 
one-body term is simply Z, the number of protons. In the RCSR, however, this term is a 
function of q and includes a kinematic factor r CT (p, q), as shown in Q4.13J ). This factor can be 
interpreted as the relativistic effect of nucleon recoil, and becomes appreciable for |q| ~M. 
In general, the evaluation of (|4.13|) requires knowledge of the full momentum distribution 
function n CT (p), which is not accurately known from experiments. This then limits the direct 
applicability of the RCSR, and in particular, the extraction of correlation information from 
the analysis. 

We therefore develop an expansion of the one-body term in moments of the nucleon 
momentum. This leads to approximate sum rules, which require only partial knowledge 
of the nucleon momentum distribution n CT (p), i.e., the first few moments, assuming that 
the moment expansion is rapidly convergent. The simplest case is RCSR-I, given in (|5.7|) , 
which involves only the zeroth moment. This version of the sum rule requires no information 
about n a (p), and therefore has the most direct correspondence with the familiar NRCSR. 
Although RCSR-I is surprisingly accurate (~ 1/10%) in approximating the one-body term 
for a uniform system of Dirac protons, it neglects the contributions from the anomalous 
magnetic moments, which turn out to be non-negligible (~23% for |q| ^> M and Z = N). 
Therefore, we keep the second moment contribution to obtain RCSR-II, given in ( |5.15| ), 
which allows the evaluation of the one-body term to within ~ 1% for a uniform Fermi gas. 
This version of the sum rule requires knowledge of only two parameters, namely (p 2 ) CT for 
protons and neutrons, which are approximately known from experimental and theoretical 
information; we have used the Fermi gas model for our estimate. This study shows that 
the uncertainty in (p 2 ) CT is likely to introduce much less error into the analysis than would 
ignoring the anomalous magnetic moments altogether. The excellent convergence of the 
moment expansion, shown here to second-order for a uniform Fermi gas, is not expected to 



be very different for correlated systems (see, for example, Donnelly et al. |16|). Furthermore, 
the extension of this approach to include higher moments, if necessary, is straightforward. 
We therfore conclude that RCSR-II, which includes contributions from anomalous magnetic 
moments, is the most efficient and reliable method for evaluating the sum rule for (e, e') data 
and separating the two-body information from the dominant one-body term. 

We now discuss the relation of our work to that of earlier authors. DeForest's|T7| approach 
to a relativistic sum rule is similar to that of RCSR-I, both formally and numerically. The 
main difference is, in the language of Section 5, his replacement of r CT (0, q) in ( |5.1| ) with the 
Lorentz invariant factor G 2 E (Q 2 ) / G 2 E (Q 2 ) = (14t)/(1+2t). However, a careful inspection of the 



recoil factor ( |4.14|) for a nucleon with arbitrary p shows that r CT (p,q) ~L o(p? q)/4E p+(l E p 
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is not a Lorentz scalar, and that anomalous magnetic moment effects must enter when Pt^O. 
Hence there is, in fact, no formal justification for writing the recoil correction at p = in 
scalar form. It is interesting that the factor (l+r)/(l+2r) is numerically nearly identical (for 
all |p| ~£>f) to our recoil factor r/(q) =r p (0, q) = (E q +M)/2E cl . However, as we have seen in 
Figure 1, the correction due to Fermi motion, which is large due to the effects of anomalous 
magnetic moments, must be included for a reliable evaluation of the one-body term. For 
example, for |q| ~ 1 GeV/c, the factor given by DeForest underestimates relativistic recoil 
corrections by ~7%, compared to RCSR-II; therefore, division of the experimental data (as 
done in Ref. |18|, for example) by DeForest's factor in the calculation of the experimental 



Coulomb sum will lead to an over enhancement of the sum by ~7%. 



The goal of Donnelly et a/.|T6| is somewhat different: to find a factor g(uo, q) which would 
replace G Ep (Q 2 ) in the definition of the spacelike Coulomb sum function, such that the sum 
will reach Z in the high-|q| limit. Although their approach can, in principle, achieve that 
goal to essentially any order of accuracy, it necessarily leads to a correction factor g(u>, q) 
which depends explicitly on the excitation energy u. However, unless the cu-dependence can 
be extracted from the Coulomb response function before integration over u, one will not 
arrive at a non- energy-weighted sum rule. In particular, the approach taken by Donnelly 
et al. will necessarily involve energy-weighted sum rules, which differ formally from the 
NRCSR, and the RCSR we have derived, in that dynamical information becomes mixed 
with the correlation information. 

A crucial step in the derivation of our sum rule is the extraction of the energy dependence 
from the response function before integration over uo. To do this, we needed two assumptions: 
first, that the off-shell form factors could be obtained from the on-shell form factors Ge{Q 2 ) 
and Gm{Q 2 ) as in ( |2.10| ). Although plausible, this requires theoretical justification, i.e., 



a theory of the electromagnetic structure of nucleons. The second assumption is that the 



form factors are all proportional, as in (|2.11|) - (|2.14|) . This is probably sufficiently accurate 



except possibly for taking Ge, u = 0, as noted earlier. A better approximation would be to 
take Ge,uO(.Ge, p , using data from larger Q 2 , although this is incorrect near Q 2 = 0. 

The results of this work do have some bearing on the question of the saturation of the 
NRCSR, which we mentioned in Section 1. The Coulomb sum function usually extracted 
from (e, e') experiments is that of ( |1.5| ), but with a finite upper limit on a>; since the theo- 
retical limit at fixed |q| is oj < |q|, this is essentially identical to ( |4.3| ). For the relativistic 
sum rules given in Section 5, one must divide the sum by the relativistic recoil factor r(q), 
i.e., either r/(q) or r//(q), to obtain a form which should "saturate" at Z in the limit that 
|q| — >oo. Since r _1 (q) > 1, this correction will enhance the experimental sum, e.g., as was 
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shown in Figures 2-4. At the energies of interest, however, this is only a small effect, e.g., 
at |q| =500 MeV/c we have r/~.94 and r//~.97. It is interesting to note that these factors 
grow with increasing (p 2 ), and therefore do show in part the increased suppression seen in 
larger nuclei. Of course, many other effects may also contribute to the observed suppression. 

Although we have worked within a "nucleons-only" approximation, in principle, N-^N 
processes are easily included in this formalism (since the matrix element of T^ between 
N states is the same as that between N states). However, we do not pursue this here 
because other N contributions, e.g., NN pair production, which may also enter the spacelike 
response, are likely to be of the same order. In order to understand fully the contribution 
of antinucleons to electron scattering, it will be necessary to study the spacelike response 
function in interacting nuclear models, and devise an appropriate renormalization scheme to 
obtain finite results. 

Two issues which are not treated in this paper must also be addressed if the relativistic 
Coulomb sum rule given here is to be used for the extraction of nucleon-nucleon correlation 
information: namely, the contribution of meson exchange currents, and internal excitations 
of the nucleon. The first can be thought of as generating two-body (or more) terms in 
the nuclear current operator J^{q) in ( |2.6| ), as done, for example in Ref. JTTJ] , where meson 



current contributions are actually calculated in a model. From the point of view of the 
sum rule, these many-body contributions can be grouped with the two-body term S^ 2 ^(q) 
of ( |4.15| ), to be determined experimentally as what is left over when £W(q) is subtracted 
from the Coulomb sum. The second issue can be thought of, to a first approximation, as the 
problem of removing the "background" of N(e, e')N* excitations from the (e, e') response, 
to leave the "nucleons-only" response, to which our relativistic Coulomb sum rules apply. 
This could be done theoretically by calculating the N* background, e.g., in a noninteracting- 
nucleus model. Alternatively, some partially exclusive experimental information could help: 
for example, (e, e'p) data in the quasi-free peak region can be used to normalize the knockout 
of protons (as opposed to N*), although other process (e.g., A— >p+7i) will also contribute to 
this process. These two problems, and the problem of antinucleons, are subjects for future 
investigation. 
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A S ingle- fermion Solutions 



The Dirac equation for free point nucleons represented by the wavefunction ip(x) can be 
written 



(z 7 ^-Af)V(x) 



0. 



(Al) 



where x M = (t, x). For a static, spatially uniform system, the solutions to (|A1|) can be 
expressed as momentum eigenstates with momentum p. The energy eigenvalues for positive 
and negative energy solutions are e p ^ = ±E p , where -E p = Vp 2 + M 2 . 

From the complete set of solutions to ([AT]), it is possible to construct a local field operator 
in terms of nucleon and antinucleon degrees of freedom: 



^(x) 



1 



E 

ps 



*MP) 



2E r 






"ps 



(A2) 



- P y — P 

where V represents the normalization volume. The nucleon and antinucleon spinors are 
defined 



««(P) 



£ P + M 






x« 



*MPJ 



£ P + M 



E p +M 
Xs 



Xs 



(A3) 



respectively, and have been normalized to 2E P particles per unit volume. Here Xs is the 
usual two-dimensional Pauli spinor, and a ps and b ps are destruction and creation operators 
for nucleon and antinucleons, respectively. 

It is possible to generalize this formalism to include isospin degrees of freedom. We 
simply define u scr (p) = tt s (p)r/ CT , where 77^ is a two dimensional spinor corresponding to 
isospin projection a, and include a sum over a in (|A2| ) . 



B Coefficients of Moment Expansion 



The expansion coefficients appearing in (|5.9|) are: 

E n + M 



'0 



2E n 



[BV 



r[ 



q 



2ME\ 



;B2) 



29 



= -q 4 (2ff q + M) (1 + ^)V (m , 

2 4M 2 £3(£ q + M) 2 2M 2 £ q (£ q + M) l ; 

P = 2E q -M(E q + M)(2ff 2 -3M 2 ) _ (1 + S ) 2 



-2 = o^wf i ^n (B5) 



4M 2 £5(£ q + M) 2M 2 £ q (£ q + M) 

n K >»q 

2 2M 2 £ q (£ q + M) 

-K 2 

^ = 2M 2 £ q (£ q + M) (B6) 

We have Tq =r™ = by our approximation G,e,„(<3 2 ) = 0. 

Expressions ( [Bo] ) and ( |BED can be obtained from QB"3| ) and (|B4[), respectively, by keeping 
only terms proportional to (1+k p ) 2 and letting (1+k p ) 2 — ^-k 2 . The corresponding expressions 
for Dirac nucleons are obtained by setting K, p = K n = 0. 
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Figure 1: Recoil correction factors vs. three-momentum transfer q. The curves representing 
rj(q) and rjj(q) with n p = K n = are nearly indistinguishable on this scale. 



Figure 2: Relativistic Coulomb sum (divided by Z) vs. three-momentum transfer q in the 
Fermi gas model, ignoring anomalous magnetic moments. The dashed curve is Si(q)/Z, and 
the solid curve is Sn(q)/Z. For comparison, the dot-dashed curve is the unmodified sum 
E(q)/Z. 



Figure 3: Relativistic Coulomb sum (divided by Z) vs. three- momentum transfer q in the 
Fermi gas model, including the proton anomalous magnetic moment. The dashed curve 
is Sj(q)/Z, and the solid curve is Sn(q)/Z. For comparison, the dot-dashed curve is the 
unmodified sum E(q)/Z. 



Figure 4: Relativistic Coulomb sum (divided by Z) vs. three-momentum transfer q in the 
Fermi gas model, including proton and neutron anomalous magnetic moments. The dashed 
curve is <S/(q)/Z, and the solid curve is Su(q)/Z. For comparison, the dot-dashed curve is 
the unmodified sum S(q)/Z. 
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